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We extend our early analyses of the dilepton and photon emission rates from a hadronic gas to 
account for strange mesons using a density expansion. The emission rates are reduced to vacuum 
correlation functions using three-flavor chiral reduction formulas, and the latters are assessed in 
terms of empirical data. Using a fire-ball, we compare our results to the low and intermediate mass 
dilepton data available from CERN. Our results suggest that a baryon free hadronic gas does not 
account for the excess of low mass dielectrons observed at CERES but do well in accounting for 
the intermediate dimuons at HELIOS. The same observations apply to the recent low and high pt 
dielectron rates from CERES. 
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OS 1 I. INTRODUCTION 

The latest experiments at the CERN SPS 0] machine have revealed a sizable enhancement of low mass dileptons 
above the two-pion threshold, triggering a number of theoretical investigations A somewhat smaller enhancement 

was also noticed in the intermediate mass region around the phi both in the dielectron (CERES) and dimuon (HELIOS) 
experiments jljjTj] . This region requires a strangeness assessment of the emission rates. An example is the recent 
analysis by Li and Gale ||. 

I , In a recent series of investigations we have assessed the dilepton and emission rates emanating from a hadronic gas 
without strangeness using general principles |],||]. We have found that the current low mass dilepton enhancement 
y—{ \ can only be accounted for by allowing for a sizable nucleon density [Q. It is important to stress that our construction 
OS ■ is not a model. It relies on the strictures of broken chiral symmetry, unitarity and data. Most model calculations 
should agree with our analysis to leading order in the pion and nucleon densities. An example being the comparison 
by Gale for the baryon free rates || . 

The aim of the present work is to extend our baryon free analysis to the strangeness sector, to account for strange 
mesons in the hadronic gas. Our results will borrow on the extension of the chiral reduction formulas to QCD with 
three flavors fl(i[| . To leading order in the meson densities, the emission rates in the hadronic gas involve forward 
scattering amplitudes of real (photon) and virtual (dilepton) photons, which behaviour is constrained by data mostly 
from electro-production and photon-fusion reactions. By including strangeness we aim at evaluating the emission 
rates in the intermediate mass region around the phi. Throughout we will not consider baryons. 

The structure of the paper is as follows : in section |J, we derive the emission rates for a baryon free hadronic 
gas to leading order in the final meson densities including strange mesons. In section |III| , we discuss the relevance 
of the various contributions and introduce the pertinent spectral functions. In section |IV| , we discuss the integrated 
dielectron and photon emission rates. In section ^, we use a fire ball scenario to account for our time-evolved rates 
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■ data. Our conclusions are summarized in section 



in comparison to current CERES and HELIOS data . In section |VT|, we discuss the current pt spectrum of CERES 
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II. DILEPTON RATES 



In a hadronic gas in thermal equilibrium, the rate R of dileptons produced in an unit four volume is directly related 
to the electromagnetic current-current correlation function pd|Jl2] |. For massive dileptons with momenta p\,2^ 
the rate per unit invariant momentum q = pi + P2 is 

where a = e 2 /An is the fine structure constant, and 

W( 9 ) = J d'xe-^Tr (e^^-F^J^O)) = i + ^ 0/T /m W F (q) . (2) 

Here eJ^ is the hadronic part of the electromagnetic current, H is the hadronic Hamiltonian, F is free energy, T is 
the temperature, and W F (q) is 
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W F (q) = i J d 4 xe iq x Ti (e-( H - F )/ T T*(J^(x)J M (0))) 
= i J dW*-»(0|T*(J'*(a;)J /1 (0)|Q) 

+ E* / / d 4 ^<7 l S,(fc)|T*(J"(x)J M (Q))K(A)) £ 



+ ••• (3) 

where the sum is over physical mesons including strange ones. For temperatures T < the first two terms in 
Eq. (||) are dominant. The first term relates to the vacuum current-current correlator and captures the essentials 
of the resonance gas model. It follows from e + e~ annihilation data. The second term is the first correction to the 
resonance gas model resulting from one meson in the final state. Two and higher meson corrections in the final state 
can be evaluated using similar arguments Q . 
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Using the definition of the electromagnetic current J^ m = q^f^Qq = + -75 V^, and the decomposition 



T*(J M J„) = T*(V^) + -T*(V»V») + -^T*(V^V« + ViVl 
we may follow P,13[ and rewrite Eq. (H) in terms of two-point correlation functions 
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ST* (V= (y)Vt(0)) 0)j = (-q 2 g^ + q^q v ) Im U c v d (q 2 ) 

5^0^(^1,(0)) 0)) =hV + ^)/mn^ 2 ), 
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and additional three- and four-point correlation functions using three-flavor chiral reduction formulas |io| . The 
axial- vector current appearing in Eq. (0) is one pion reduced p,[l^]. 
With the above in mind, we have 



ImW F (q) = -3q 2 Im ^(q 2 ) + + J j^W[(q,k), 



(6) 



with Hy = Ily 3 and Hy = flly 3 . The first term is the analogue of the resonance gas contribution with no chiral 
reduction involved. The second term is the correction to the resonance gas model resulting from one meson in the 
final state. Use of the on-shcll three flavour chiral reduction formulas gives for the part involving solely the two-point 
correlators 
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Wf (q, k) = gl-^q'lm Il v (q 2 ) + g^q 2 Im U v (q 2 ) + -IL v (q 2 ) 
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-9l T2 {k + qflm n z A ({k + q) 2 ) ~ g« '73- (* + qf (im I%{{k + qf) + Im Tl U A ((k + q) 2 )) 

J it J, 



K 



+9k^ - ™W)Re (A£(* + q) + A£(fc - q)\ Im n v (q 2 ) 
+ 9 £j r (v 2 -mW)Re (AK(k + q)+A%(k-q))lm (u v (q 2 ) + h%(q 



where v = k ■ q. The contributions from three- and four-point correlators read 



Wf( 9 ,A) = ^jf3V3- 
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+permutations(g — > —q, k — > —k) 
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where the permutation applies only for the last term, and no-mixing is assumed ( n 3 / = Ily = ). The thermal 
meson density function g k and the meson propagator A^(g) are defined by 
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The three point correlator is given as 



z Jy 



ST* 



V^ 3 (z) + -±=V^ s (z)) (v 3 Jy) + iv»l ^(0) 



V3 



V3 
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The meson decay constants, Jk ~ 1.24/ T , f v w 1.32/^, /,,/ w 0.74/ ff are used ||14|. The results Eqs. (|6f||) reduce to 
the two-flavor results discussed in Since they include strangeness they can be used all the way through the phi 
region. 

The spectral functions appearing in Eqs. (^,0) are related to e + e~ annihilation data. They will be borrowed from 
experiment. The three- and four-point correlation functions VVj^j^, Vj^j^ and VVcr are constrained by the two- 
photon fusion reactions and crossing symmetry [fl0| . A detailed analysis of these processes show that their contribution 
to our rates can be ignored for both CERES and HELIOS (2% correction). These observations confirm and extend 
the ones made in || to the three-flavour case. Hence, the two- and three-point contributions will be omitted for most 
of the discussion to follow. 



III. SPECTRAL FUNCTIONS 



In this section we derive explicit expressions for the spectral functions. The empirical information on the two-point 
correlators will be inserted in the rate production by means of suitable spectral weights. For instance, the contribution 
of the non-strange hadrons to the spectral weights is dominated by the p (m p = 768.5 MeV, To.p = 150.7 MeV) and 
the ai (m ai — 1230 MeV, To tai = 400 MeV) resonances. For that, we will follow the arguments of || and use the 
following parametrizations 



fp m l + 7g 2 

q 2 m 2 — q 2 — im p T p {q 2 ) 



m«! - Q 2 - im ai T ai (q 2 ) 



with f p = \[2f v and f ai = 190 MeV. The decay widths T are given by 

,2 _ Arm 2 \ 3/2 
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i 2 p ~ Am 2 
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where 9(x) is Heaviside functions (step function). Analogous parametrizations will be used for the strange resonances 
as well including the w and <j> in Hy, and K\ in II4. The isovector part of the spectral weight II y is dominated by 
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the p and its radial excitations, while the hypercharge part Tly is dominated by the omega, the phi and their radial 
excitations. In the physical basis 




nf? = in£ + |n$ . (13) 

In Fig. |l| we show the electromagnetic spectral function Ily™ = Hy + j^ly following from our parameters (solid 
curve) in comparison with the data compiled in [flEf . 
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FIG. 1. Electromagnetic Spectral Function. See text 

The resonance parameters we have used are summarized in Table. || |l6| . The decay constants are fit to the empirical 
spectral weights [JL4|. They are all within 10 % of the constituent quark model. Although the axial- vector spectral 
weight is not well-known at large invariant mass, we observe that its numerical contribution to the dilcpton emission 
rates is overall small. 







I G (J PC ) 


mass (m,i) 


decay width (G;) 


decay constant (fi) 




P (770) 

p(1450) 

,9(1700) 


!+(!-) 


768.5 
1465 
1700 


150.7 
310 

235 


130.67 
106.69 
75.44 




w(782) 

w(1420) 

w(1600) 


o-(i-) 


781.94 
1419 
1649 


8.43 
174 
220 


46 
46 
46 




0(1020) 
0(1680) 


o-(i-) 


1020 
1680 


4.43 
150 


79 
79 


ni 


oi(1260) 


1-(1++) 


1230 


400 


190 (/„) 




A"i(1270) 
Ki(1400) 


^(1 + ) 


1273 
1402 


90 
174 


90 
90 



TABLE I. Resonance parameters 
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The high energy tail of the vector spectral weight is still above the free qq threshold well through the high mass 
region as shown in Fig. for i?>3 = 127rlmil(s) without the J '/SP contribution. We will use the parametrization 



0.9 , 

— 1 + tanh 

8x \ 



(14) 



as shown by the solid line. Here q\ = 2 GeV and q = 0.4 GeV. We recall that the free qq threshold for SU(3) is 



*>a = £ 



<Ti(> 3tt) 



a(e + e — > fi+fj, ) 



2 ^ ~ q 

q—u : d,s 



= 2 



(15) 



and underestimates the hadronic correlations by about 30% in the 2.5-3.5 GeV region. This tail is important in the 
emission rates as we will note in the next section. It maybe interpreted as (logarithmic) corrections to the perturbative 
quark result by duality. Above 3.5 GeV the charm effects show up. 
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FIG. 2. R >3 values. Data from Adone 17fl and SLAC-LBL tg 



IV. EMISSION RATES 



A. Dielectron Rates 



Given the two-point spectral weights, it is then straightforward to reconstruct the emission rates using Eqs. 
The contributions Eq. (j8|) are found to be numerically small and will be ignored. R can be reexpressed in terms of 
the invariant dielectron mass M = yfq*, the rapidity rj and the magnitude of the transverse momentum q t fll|| 



dR 

d 2 M 



= / dy 



x dR 

2 



dq 2 t t: ^{M,ri,qt). 



(16) 



In Fig. ||, we show our new three flavor-rate (solid line) in comparison to the two-flavor rate || (dashed line) . Aside 
from the omega and phi which were absent in the two-flavor analysis, the results are overall consistent for T = 150 
MeV. The chiral reduction results show a substantial enhancement of the low mass dielectrons in comparison to a 
simple PC AC treatment [§. 
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FIG. 3. Dielectron rates at T = 150. See text. 



In Fig. ^ we show the rates from a hadronic gas in comparison to a free quark rate from a quark gas C^] at T = 150 
MeV. The plasma rates are large in the low mass region (more than a factor of 2), above the phi (more than a factor 
of 2), and around the 2 GeV region (about a factor 1/2). Surprisingly, however, the hadronic tail Eq. ( O ) (e + e~ into 
hadrons) still provide substantial emission strength even through the high mass region in comparison to a free quark 
gas at the same temperature. An enlargement of that region is shown in Fig. || 
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FIG. 4. Hadronic gas (bold) versus quark gas (thin). See text. 
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FIG. 5. Same as Fig. §. 



B. Photon Rates 

The rate for photons follows from Eq. (|^) at q 2 = 0, 

<°lr-^ w <"- (17) 

Fig. [| shows the photon emission rate for T — 50, 100, 150 MeV. They are overall consistent with the photon rates 
established for the SU(2) case ||. 
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FIG. 6. The photon emission rate from a hadronic gas. 
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V. EMISSION RATES FROM A FIRE BALL 



In a heavy-ion collision the electromagnetic emission occurs from various stages of the collision process. In this part, 
we will focus on the low and intermediate dilepton emission rates (up to 1.5 GeV). We will assume that they emanate 
from a simple fire ball composed of a hadronic gas. The fire ball will be modeled after transport codes B, for the 
CERES (S-Au) and HELIOS-3 (S-W) heavy-ion collisions. In particular, the expansion will be assumed homogeneous, 
with a volume and temperature parametrized as J3],[3| 

V(t) = Vo ( 1 + * 



to, 

T(t) = (T t -T 00 )e- t ' T +T X (18) 

with t = 10 fm/c, T, t = 170 MeV, = 110 MeV, t = 8 fm/c, and the value of V is absorbed into the over-all 
normalization constant N Vq (= 6.76 x 10 -7 fm 3 ) in Eq. (|l9|), which is fixed by the transport results. The freeze-out 
time will be set to i/. , = 10 fm/c. Since we would like to only address the issue of a thermal gas of mesons, we will 
ignore the baryons assuming their contribution at present energies to be small in the leading density approximation [0 . 
The mesonic fire ball is more appropriate for RHIC as opposed to CERN energies. 



Using Eq. (18) and the above rate, the final expression for the integrated emission rate per unit rapidity r) and 



invariant mass M is 



dN/dndM,^, Ar , r [**•'■ , T „, fd 3 q A , 2 ,dR 

-(M) = N M dtV(t) —±A(q ,q 2 )— . (19) 



dN c h/dr) J J q a d 4 q 

The acceptance function A(q ,q 2 ) enforces the detector cut p t > 200 MeV, 2.1 < rj < 2.65, and 6 ee > 35 mrad 
for CERES. For HELIOS-3, Eq. ( |i~9| ) can be used by integrating over rj with the cut m t > 4(7 — 2rf),m t > 
v / (2m M ) 2 + (15/cosh(r/)) 2 . The results for CERES and HELIOS-3 are summarized in Figs. @ and § 
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FIG. 7. Dielectron rates for CERES S-Au experiment. See text. 
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FIG. 8. Dimuon rates for HELIOS-3 S-W experiment. See text. 

The thin line in Fig. ^ is our result using only the two-point spectral weights in Eq. , while the thick line follows 
from an additional gaussian smearing over the detector mass resolution, 

' ,„r/ 1 ( i M ~ M') 2 \ dN/dridM ^ r/ . 



with a — 0.1 x M. The background (dotted line) was taken from a transport model [g|JS|]. The 4> resonance is clearly 
visible. The data are from [1|. Clearly the fire-ball evolution together with the mesonic gas do not account for the 
low mass dileptons. In Ref. [3], it was shown that only large nucleonic densities could account for the data to leading 
order in the hadronic densities. 

The solid line in Fig. [8] is our three-flavor result for the HELIOS-3 experiment. The dashed curve refers to the 
results discussed recently by Gale and Li || using effective Lagrangians and a variety of two-body reactions. Our 
results are in overall agreement with theirs. The dotted line is again the background contribution from the transport 
model The fair agreement of the fire ball with the dimuon spectra above 1 GeV indicates that a thermal 

hadronic gas treatment is consistent with the data. Some enhancement in the low mass region maybe achieved by 
adding baryons, but not enough in our leading density approximation to bring it into agreement with the data pfl. 



VI. CERES WITH P T CUT 



Recently, the CERES collaboration analyzed the pt dependence of the dielectron pairs [g0| . Using our mesonic rates 
and the fire ball evolution Eqs. (jig- 19) we show in Figs. [)] and |lC] our results in comparison to data. The background 



contribution in our case was borrowed from transport calculations pl[ . The initial temperature of the fire ball is set 
at Ti — 160 MeV and the freeze out temperature at = 105 MeV for which the time scale is tf. a , = 20 fm/c, with 
r = 10 fm/ c and to = 10.8 fm/c. The new normalization constant, NoVq = 3.45 x 10~ 7 fm 3 , is fixed by the transport 
results |lj. 

Our analysis based on a baryon free hadronic gas (with strange mesons) does not reproduce the low mass dielectron 
enhancement in both Fig. |^ and [l^. Since most of the discrepancy of the low mass enhancement comes from low p t 
contribution, taking into account the large statistical and systematical errors for M > 1.0 GeV, we find that our high 
Pt spectrum is consistent with the data. In Fig. [ll], the data are the differences between the mean values of the total 
and low p t spectrum, so we do not include the vertical error bars. 
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FIG. 9. Dielectron rates for the total pt of CERES Pb-Au experiment. 




FIG. 10. Dielectron rates for low p t of CERES Pb-Au experiment. 
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FIG. 11. 



Dielectron rates for high pt of CERES Pb-Au experiment. 



VII. CONCLUSIONS 



We have extended the dilepton and photon emission rates from a thermalized gas of hadrons to account for strange 
mesons without baryons. The analysis relies on the extension of the chiral reduction formulas from SU(2) to SU(3). 
The hadronic gas now includes the effects from kaons, etas and phis. As expected, only substantial changes in the 
rates are seen around the omega and phi region. 

In the low mass region there is a substantial enhancement due to the mixing with the axial particles, but not enough 
to account for the CERES data after time-evolution and detector cuts. In the intermediate mass region, there is good 
agreement with the HELIOS data. The baryonic effects are important in the low mass region, but not enough to 
account for the discrepancy seen in the Pb-Au collisions at CERN within the density expansion ||. The large and 
persistent hadronic tails make the thermal mesonic emission rates comparable to the ones from a thermal gas through 
the intermediate mass region. 

As we have indicated in the introduction, our analysis is not based on a model. Our only assumption is that the 
emission rates follow from a baryon free and dilute hadronic gas. The rest follows from symmetry and data. As such 
all model calculations should be in agreement with our results under these assumptions, and it is gratifying to see 
that the recent analysis performed by Gale and Li || using reaction rates does. 

The fact that the HELIOS-3 data can be explained without recourse to additional assumptions, make part of our 
arguments reliable. The persistent disparity in the low mass region with CERES and the newly measured pt spectra 
may indicate the need for an enlargement of the original assumptions. The importance of the thermal quark gas 
emissivities in the 1/2 and 1 GeV region maybe indicative of some simple n onp erturbative effects in the partonic 
phase for the temperatures considered (22). This issue will be addressed next ]l9|. 
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